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Study of Mrk 501 above 60 GeV with CELESTE 

E. Brion, for the CELESTE collaboration 

CENBG, Chemin du Solarium, Le Haut-Vigneau, BP 120, 33175 Gradignan Cedex, France 

The CELESTE atmospheric Cherenkov detector, running until June 2004 at the Themis solar 
facility, has taken data on compact sources such as pulsars and blazars. We will take stock of the 
experiment, in particular regarding the latest improvements of the detector simulation and data 
analysis. These changes provide us with a new analysis of old data with smaller uncertainties. We 
present here the evidence for a weak signal from Mrk 501 in 2000-2001. 
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I. INTRODUCTION 

CELESTE (Cherenkov low energy sampling and 
timing experiment) was a Cherenkov experiment u- 
sing the heliostats of the former Electricite de France 
solar plant in the French Pyrenees at the Themis site. 
It detected Cherenkov light from electromagnetic sho- 
wers produced in the atmosphere by the 7-rays coming 
from high energy astrophysical sources. The light is 
reflected to secondary optics and photomultipliers in- 
stalled at the top of the tower. Finally it is sampled 
to be analysed [jj. 

Two states of the experiment have to be distin- 
guished to classify CELESTE data. During the first 
one (between September 1999 and June 2001), 40 he- 
liostats were used with two types of pointing (single 
pointing: all heliostats at 11 km, double pointing: 
half at 11 km, half at 25 km), and during the second 
one (between September 2001 and June 2004), 53 he- 
liostats pointing at 11 km were used (of which 12 veto 
heliostats aiming wide for proton rejection). An ana- 
lysis improvement has been made on the last data with 
53 heliostats taking the Crab (Nebula) data as refe- 
rence, since the source is bright and stable. The new 
analysis variable provided a sensitivity of 5.1 cr/Vh. on 

5.3 h data, whereas the old data analysis gave 2.0 and 

3.4 a/Vh for single and double pointing Q. The new 
analysis has been tested on old Crab data with 40 he- 
liostats and gives better sensitivity as will be shown. 

CELESTE has taken 40 heliostats data on the 
blazar Mrk 501. The old analysis gave 2.5 a We 
present here the results obtained with the new ana- 
lysis. Mrk 501 is a source interest: it is one of the 
brightest and closest (z = 0.034) extragalactic X-ray 
and TeV sources. A detection with CELESTE could 
help to understand the emission processes. 

First, the new analysis will be explained with the 
results obtained on the Crab data. Then, the data 
selection and analysis, and finally the light curve of 
Mrk 501 will be presented. 



detect low energy 7-rays, we use the sum of the in- 
dividual digitized signals to increase the signal-to- 
background ratio (figure nj. The summation includes 
a correction for the sphericity of the wavefront, as- 
sumed to be centered in the 11 km plane. Assuming 
a wrong position for this center (impact parameter) 
broadens the sum: the height-over- width ratio (H/W) 
decreases. We compute (H/W) repeatedly for a grid 
of different assumed positions. The impact parameter 
is taken to be the position for which the (H/W) ra- 
tio is maximum, denoted by (H/W) max . This is valid 
for 7-rays (figure (a)) but not for protons arriving 
isotropically and for which the wavefront is not sphe- 
rical (figure El (b)). A measurement of the flatness of 
these 2D-distributions is given by the following esti- 
mator: 



aV6rage { (H/W)^ 



over 24 positions 



where (H/W)2oom is the average of (H/W) over 24 
positions along a ring 200 m from the maximum po- 
sition. For 7-rays there is a clear maximum and this 
estimator takes low values. For hadrons, it is usually 
larger (figure (c)). 

To improve the signal quality the cut was optimized 
on the Crab data. For the data with 53 heliostats, the 
applied cut is £ < 0.35 (figure EJ) which is the biggest 
contribution to an improvement in Crab sensitivity 
to 5.1 cr/Vh. The detection is stable. 40 heliostat 
Crab data were also analysed with this variable. Since 
this estimator is very sensitive, the chosen cut is diffe- 
rent: optimizing the cut on 40 heliostats with single 
pointing gives £ < 0.30 (£ < 0.35 for double poin- 
ting) which is not yet understood. Simulations are in 
preparation to confirm (or not) this cut. The results 
are presented in table [I] 



III. DATA 



II. ANALYSIS IMPROVEMENT 

A new method was found to exploit the FADC in- 
formation to reject the hadronic background pj. To 



CELESTE took data on Mrk 501 during 2000 and 
2001 with 40 heliostats. The previous analysis giving 
2.5 a, we hope to improve this result with the new 
analysis. 
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Y: spherical wave front 
- centered in 

the 11km plane 




high & thin 
FADC sum 



low & wide 
FADC sum 



FIG. 1: The shape of the FADC sum depends on the po- 
sition of the center of the wavefront we assume |5|. 



(a) Gamma 100 GeV 



(b) Proton 500 GeV 






FIG. 2: Height-over- width ratio (H/W) computed for dis- 
crete positions in the 11 km plane (a) for a 100 GeV 7-ray 
and (b) for a 500 GeV proton, (c) Distribution of the £ 
parameter for OFF data (filled black circles) and 7-rays 
simulation (red line). 



A. Data selection 

Data selection is stricter than in We keep 

data with stable currents, small ON— OFF current 
and trigger rate differences and number of events 
> 10, 000. Finally, after analysis cuts, the OFF data 
rate should only depend on the hour angle. Low rates 
at small hour angle are an indication for bad atmo- 




FIG. 3: Significance a on the 53 heliostat Crab data in- 
tegrating f by the left: 1(f) = J crdf. The maximum of 
the signal is obtained for £ < 0.35. 



Data set 


Applied 


Significance 


Duration 


Sensitivity 


Flux 




cuts 


a 




[a/Vh] 


[7/min] 


40 heliostats 


£ < 0.30 


13.0 


8.4 


4.5 


2.9 ±0.2 


single pointing 












40 heliostats 


f < 0.35 


11.9 


8.4 


4.1 


5.3 ±0.5 


single pointing 












40 heliostats 


£ < 0.35 


9.7 


5.7 


4.1 


5.0 ±0.5 


double pointing 












53 heliostats 


£ < 0.35 


12.5 


5.3 


5.4 


4.5 ±0.4 


with veto 












53 heliostats 


f < 0.35 


11.6 


5.3 


5.1 


3.5 ±0.3 


with veto 


no veto 











TABLE I: Sensitivity and flux corrected from hour angle 
efficiency for the different data sets on Crab. All data have 
an absolute hour angle less than 1.5 h. 



spheric conditions. Thus, we require that they be 
> —7.5 x |ah| + 15 (figure^). There was no cut applied 
on hour angle. 



B. Data set and analysis 

After selection, 37 pairs (10.9 h) remain: 10 pairs 
were taken with single pointing at 11 km and 27 with 
double pointing at 11/25 km. All pairs are kept to- 
gether to have more statistics. In terms of years, we 
have 31 pairs during 2000 and 6 pairs during 2001. 

The results of the new analysis are presented in ta- 
ble [H] Since 40 heliostat Crab data with single and 
double pointing don't give the same cut to apply on 
£, the difficulty here is to choose the cut (figure [5} . If 
we apply £ < 0.30, we get 2.9 a on all the data (3.1 a 
in 2000, 0.0 a in 2001). If we apply f < 0.35, we get 
again 2.5 a on all the data like with the old analysis 
(2.4 a in 2000, 0.7 a in 2001). 

The chosen cut to present the lightcurve hereafter 
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1.5 2 2.5 

Absolute value of hour angle [h] 



FIG. 4: Mrk 501 OFF data rate after analysis cuts as a 
function of hour angle (typical raw trigger rate is 20 Hz) . 
Data are taken with 40 heliostats single pointing (green 
circles) and double pointing (red triangles). 
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FIG. 5: ON-OFF f distribution for the 40 heliostat 
Mrk 501 data (single and double pointing). 



is £ < 0.30, because we have the assumption that the 
distribution of the £ variable depends on the declina- 
tion (owing to its construction). Figure shows this 
distribution for the 53 heliostat Crab data. In com- 
parison, we can see in figure [7| the same distribution 
for 53 heliostat Mrk 421 data, which has the same 
declination as Mrk 501. The distribution is shifted 
towards the left. 



Data set 


Applied cuts 


ON 


OFF 


ON-OFF 


Significance 


Signal-to-noise 












a 


ratio [%] 


All data 


raw trigger 


789,648 


787,452 


2,195 


1.5 


0.3 




£ < 0.30 


10,556 


10,070 


486 


2.9 


4.8 




f < 0.35 


44,218 


43,363 


856 


2.5 


2.0 


Year 2000 


raw trigger 


666,704 


663,474 


3,230 


2.3 


0.5 




£ < 0.30 


9,424 


8,940 


485 


3.1 


5.4 




£ < 0.35 


38,753 


37,977 


776 


2.4 


2.0 


Year 2001 


raw trigger 


122,944 


123,978 


-1,035 


-1.9 






f < 0.30 


1,132 


1,130 


1.8 


0.0 


0.2 




£ < 0.35 


5,465 


5,385 


80 


0.7 


1.5 



TABLE II: Mrk 501 number of events (40 heliostats single 
and double pointing) after different cuts. Note that the 
number of events is rounded. 



A work is in preparation that compares thoses dis- 
tributions for Crab, Mrk 421 and Mrk 501 40 heliostat 
data with simulations. Since we have a lot of data on 
Mrk 421, we will be able to optimize our cut on this 
source and apply it to Mrk 501. 



C. Seasonal changes 

CELESTE raw trigger rates are typically 25 Hz du- 
ring the Winter. In early Fall and in Spring, rates fall 
by a factor of two. The LIDAR on site showed that 
increased atmospheric extinction due to aerosols is a 
major cause Q- Unfortunately, it started to get data 
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FIG. 6: ON-OFF f distribution for the 53 heliostat Crab 
data (filled black circles) compared to simulations on Crab 
normalized to the data (red line). 



in December 2001 after these data on Mrk 501. We 
looked at the possible influence of this effect on the 
ON— OFF excess, but nothing was highlighted: for 
trigger rate < 15 Hz, 20 pairs give 2.7 a during 6.0 h, 
and for trigger rate > 15 Hz, 17 pairs give 1.5 a during 
4.9 h. 



IV. LIGHTCURVE 

We present in figure |S| the lightcurve obtained by 
CELESTE for monthly averages with the chosen cut 
£ < 0.30. Note that Mrk 501 was active in X-rays du- 
ring 2000 and quiet during 2001 pi|. During the year 
2000, the cut to be confirmed, we get an ON— OFF 
excess of 3.1 a during 9.1 h. 
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CONCLUSION 
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1500 
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Before the LIDAR study, the acceptance was un- 
known for data with trigger rates less than 15 Hz. We 
now have more confidence in the 3.1 a obtained for 
the year 2000 and in our acceptance. Therefore, if 
this result is confirmed, we will interpret this excess 
as a weak signal and calculate a flux. The acceptance 
will include the degraded atmospheric extinction. 



FIG. 7: ON-OFF f distribution for the 53 heliostat 
Mrk 421 data 0. 
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If we interpret this excess as a signal, and under 
the (false) assumption that we have the same detector 
acceptance and spectrum for Mrk 501 as for the Crab, 
we get one fourth of the Crab flux. This is comparable 
with the results presented in 9] (data on April 29, 
1998 in X-rays similar to Spring 2000). 



FIG. 8: Mrk 501 lightcurves seen by CELESTE above 
60 GeV in 2000-2001, RXTE/ASM in X-rays 8], and 
the optical telescope of Bordeaux 0. The CELESTE 
light curve is uncorrected for efficiency in hour angle, nor 
for different atmospheric extinctions. 
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